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Abstract 

Meshes, as used in the finite element method, can be viewed as shght generahzations 
of polytopal complexes, and the volume function for tetrahedra can be generalized to a 
volume measure for polyhedra and, more generally, polyhedral meshes. Here, polyhedral 
meshes are considered up to translation and scaling, and the resulting quotient space is 
diffeomorphic to a sphere. If we choose an embedding of this sphere in the Euclidean space of 
all polyhedral meshes, then the volume function restricted to the sphere is a quality measure 
for these polyhedral meshes. In the finite element method, one is usually concerned with 
tetrahedra, hexahedra, prisms and pyramids as building blocks for meshes. We show that the 
(normalized) gradient flow of the volume is a fast and efficient regularization scheme known 
as the geometric element transformation method (GETMe). As part of this discussion, we 
prove that the volume is a Morse-Bott function for tetrahedra. 

1 Introduction 

Meshes, as used in the finite element method, can be viewed as slight generalizations of poly- 
topal complexes. In this context, mesh quality affects numerical stability as well as solution 
accuracy of this method [18]. The geometric element transformation method (GETMe) as 
a regularizing method for tetrahedra and a smoothing algorithm for tetrahedral meshes is 
introduced in [22] and generalized to other polyhedral meshes in [20]. Several numerical 
tests have confirmed that the geometric element transformation reliably and efficiently reg- 
ularizes the polyhedron types, which are relevant for the finite element method. It also 
significantly reduces errors in solutions of the finite element method and improves solution 
efficiency [21]. However, there was no mathematical proof, why GETMe works. Providing 
such a proof is the main motivation for this work. We develop the mathematical foundation 
for it and thereby discover a beautiful instance of Morse-Bott theory by way of the (signed) 
volume function. We can generalize the volume function to a (mean) volume function for 
all polyhedra and polyhedral meshes. This is clearly a scalar function of node positions 
that measures some geometric property of the element, which Knupp [11] defines to be a 
quality metric. Let us reserve the term metric for the Riemannian metric introduced on the 
manifold of all polyhedra in Section 2, and instead call this kind of scalar function a quality 
measure. For a collection and mathematical discussion of quality measures we refer to [11] 
and the references therein. 

At first sight, the volume seems like a rather useless quality measure, because it lacks 
scaling-invariance and cannot capture information about the (scaling-invariant) shape of the 
object. It is therefore surprising, that there should be a simple way to turn it into a useful 
quality measure. To this end, we introduce a novel approach to constructing and describing 
quality measures by 



1. modeling the space of polyhedra and polyhedral meshes (up to translation and scaling) 
as a differentiable manifold diffeomorphic to the sphere, 

2. embedding this sphere in the Euclidean space of all polyhedra and polyhedral meshes, 
and 

3. considering the restriction of the volume (or some other geometric function) on poly- 
hedra and polyhedral meshes to the sphere as a quality measure. 

Note that this approach easily generalizes to polytopes and polytopal meshes in arbitrary 
dimensions. The following main properties differentiate the resulting quality measures from 
other quality measures [11]: 

• The dependence of the measure on the polyhedron as a whole is obvious. This allows 
for a systematic development of other quality measures, in particular problem specific 
ones, by changing the underlying geometric function or the embedding of the sphere. 

• The measure is not only defined for valid meshes, but for all meshes. The application 
of the gradient flow, that is, the flow of the gradient vector field, of the measure 
restricted to the sphere will — by the very nature of the gradient — improve the quality 
of the polyhedron or the polyhedral mesh. The discretization of the flow yields a 
transformation method, which is expected to converge to (local) maxima of the quality 
measure. 

We work out the measure and its gradient in one specific example: With the choice of em- 
bedding from Section 2, the volume function turns out to be a simple regularity and quality 
measure for the main building blocks for meshes, namely tetrahedra, hexahedra, prisms, 
pyramids. Furthermore, it yields and generalizes GETMe [22] as a smoothing algorithm for 
meshes. Besides the fundamental importance of our results for applicational aspects of the 
finite element method, the techniques demonstrate a fruitful direction of research using the 
interplay of global analysis, discrete geometry and dynamical systems. 

In the context of tetrahedral meshes [22] , the geometric element transformation is based 
on shifting vertices by the opposing face normals. We will see in this work, that, up to 
translation and scaling, the space of tetrahedra can be viewed as an 8-dimensional sphere, 
the tuples of face normals as the gradient field of the volume function on this sphere and 
the smoothing as the gradient flow. It turns out, that this flow has only one sink, namely 
the positively oriented regular tetrahedron. The volume extends to a volume on tetrahedral 
meshes by considering the sum of the volumes for all tetrahedra in the mesh. This idea 
generalizes to a mean volume on a sphere of polyhedra and polyhedral meshes by averaging 
over all possible triangulations of the polyhedra. In particular, this approach shows the 
regularization property of the transformation corresponding to the gradient held for certain 
polyhedra and polyhedral meshes. 

In Section 2 the structure of polytopal complexes and meshes is described. Section 3 in- 
troduces a convenient notation for face normals and uses it to describe the geometric element 
transformation for a single tetrahedron. Section 4 contains a qualitative discussion of certain 
gradient flelds as well as their normalization and shows that the face normals of tetrahedra 
yield the gradient field of the volume function. We generalize these gradient fields to all 
polyhedral meshes in Section 5 and compute them for a few fundamental polyhedra. This 
naturally leads to the notion of symmetry in Section 6 for polyhedra and a new regularity 
and quality measure for polyhedra and polyhedral meshes given by a mean volume function 
in Section 7. In Section 8 we show that the volume function for tetrahedra is a Morse-Bott 
function, and Section 9 discusses the dynamics of the gradient flow for general polyhedra. 
We summarize our results in Section 10. 

2 The manifold of polyhedra and meshes 

Deflnitions for polytopes and meshes vary in the literature. For our purpose, one is as good 
as any. We base our deflnitions on [24]. 

A poly tope P is the convex hull of a flnite set of points in some R''. A linear inequality 
cx < Cq is valid for P, if it is satisfied for all points a; G P. A face is any set of the form 

F = Pn{xeIV^\cx^ Co}, 
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98 where cx < co is a vahd inequahty for P. The dimension of a face is the dimension of its 

99 affine hull. The faces of dimension 0,1 and rf — 1 are also called vertices, edges and facets, 
respectively. An isomorphism between two polytopes P C R'' and Q C R"^ is an affine 

101 map /: R"* — > R*^, that is a bijection between the points of the two polytopes. Note, that 

102 / : R'' — R*^ itself does not need to be injective or surjective. A polytopal complex C in R'' 

103 is a set of polytopes in R'' such that 

104 1. the empty set is in C, 

105 2. for any P G C all faces of P are in C, 

106 3. the intersection of any two polytopes in C is a face of both. 

107 A complex I? is a subcomplex if I? C C Let \C\ := Upec maps between polytopal 

108 complexes are all the maps /: \C\ — ^ \'D\, whose restriction to polytopes in C are afhne 

109 linear, up to isomorphism for the polytopes of C and T). Two polytopal complexes C and 2? 

110 are isomorphic, if there is a map f : C ^ V which is a bijection between C and V — {f{F) \ 

111 F G C}. The complex C{P) of a polytope P is the complex of all faces of P. A mesh J\4 is 

112 given by the pair {C,p), where p: V{C) — > R'^ is a map from the vertices V{C) of C to R''. 

113 A mesh (2?, q) is a submesh of (C,p), if I? is a subcomplex of C and q = p\v{c)- Two Meshes 

114 iC,p) and {'D,q) are isomorphic, if C and 2? are isomorphic. 

115 In this work, we will focus on three-dimensional polytopes, which — in contrast to [24] — 

116 we call polyhedra. The graph of a polyhedron is the graph of its vertices and edges. By 

117 Steinitz's theorem [17, 9, 24], the polyhedral graphs are exactly the 3- vertex-connected 

118 planar graphs. A graph is planar, if it can be embedded in the plane, and it is 3-vertex- 

119 connected, if the graph remains connected, when you delete fewer than 3 vertices from the 

120 graph. From now on, all complexes C and meshes (C,p) are polyhedral, that is, C — [J^C{P^) 

121 for some polyhedra P* and p: V{C) ^ It^, and we call a mesh (C,p) a polyhedron mesh, if 

122 C = C{P) for a polyhedron P. 

123 We are usually interested in transformations of a mesh that do not change the isomor- 

124 phism class. If the isomorphism class is fixed, we will fix a polyhedral complex C in this class 

125 (and therefore a choice of orientation) and an ordering of vertices V{C) — {Vi, . . . , Vn}, which 

126 allows us to identify the set of meshes {C,p) in this isomorphism class with points in Eu- 

127 clidean space R^" = (R'^)" via (C,p) >-)■ {p{Vi), . . . ,p(y„)). In this case let us abuse the no- 

128 tation and denote the mesh by the sequence of vectors p = {pi, . . . ,pn) = {p{Vi), . . . ,p{Vn)) 

129 and all polyhedral submeshes by (necessarily orientation preserving) subsequences of p. 

130 Notice that Steinitz's theorem allows us to alternatively define a mesh to be a 3-vertex- 

131 connected planar graph together with an embedding of its vertices into R"^. 

132 Fix an isomorphism class of a mesh. We want to consider the meshes up to translation 

133 and scaling. More precisely, consider M — R'^"\{(po, . . . ,po) \ pq S R"^} and its quotient 



N = M/ r^, where 




for A > and 
,Po)+P forpoeR^. 



Let us identify N with S^^^ i) i = 5'3(" i) i x {0} C M via the homeomorphism induced 
by the smooth map tt = a o t: M S*"^*^"^^-*^^, where 

r: R3" -> R3("-1) X {0}, {pi, . . . ,Pn) ^ {pi - Pn, ■ ■ ■ ,Pn-l - Pn,0) 

R3"\{0}^53„-i^ P^P. 

\\P\\ 

134 Consider the differential structure and Riemannian metric on N pulled back via tt from the 

135 standard 5'3("-i)-i. Notice that Tr{p) = r(p) = cr(p) = p ior p € N. We canonically identify 

136 the tangent space TpM of M at p with (R'')" and denote the R^-component functions of a 

137 (tangent) vector field X on M by ATi, . . . , Xn '■ M — > R"^. Vector fields on M/ ^ are of the 

138 form D{Tr){X) for some vector field X on M. For notational convenience, we will use column 

139 or row vectors interchangeably. Despite its simplicity we have not found this useful model 

140 for meshes anywhere in the literature. However, it is essentially the same as the model for 

141 polymer used in [4]. 
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3 Face normals 



143 Recall that the geometric element transformation for tetrahedral meshes was constructed 

144 in [22] by shifting the vertices of each tetrahedron by the opposing face normal, normalized 

145 so that the transformation is scaling invariant, and then averaging the shifts over the whole 

146 mesh. Let us introduce a convenient notation, which helps us to organize the information 

147 about face normals for general polyhedra and meshes. For (pi, . . . ,Pk) C (R"^)*^ we define the 

148 closed 1-chain C(l, . . . , fc) as the sum of 1-simplices between pj and Pj+i for j = 1, . . . , fc— 1 

149 as well as between pk and pi, as indicated in Figure 1. Furthermore we let 



. . . X p2 +P2 X p3 + . . . +pk-i y-Pk+Pk X _pi. 




Figure 1: The chain C(l, . . . , fc) 



150 There are several useful facts about cross products, which can be proven by elementary 

151 calculations. 

Lemma 3.1. We have the identities 

(3.1) X (P2-P3) = K1,2,3), 

(3.2) (pi - P3) X {p2 - pa) = 2, 3, 4), 

k-l 

(3.3) v{2, . . . , fc, 1) - ...,k) = J2 H^J^ J' + !)• 

152 The direction of the normal vector j/(l, 2, 3) can be found by a right-hand (grip) rule. Equa- 

153 tion (3.3) encodes the following properties of cross products: 

154 1. The face normals of tetrahedra add up to zero. 

155 2. If C(l, . . . , fc) is the boundary of a sum of 2-simplices, then the sum of their normals is 

156 independent of the choice of 2-simplices. Here, the sum is a signed sum, where edges 

157 with opposite signs cancel. 

158 3. If C(l, . . . , fc) is planar and convex, then \\v{l, . . . , fc)|| equals twice the area enclosed 

159 by C(l, . . . ,k). 

160 The above lemma gives a convenient way of organizing linear combinations of face nor- 

161 mals. In particular, it allows us to write the (unnormalized) geometric element transforma- 

162 tion applied as a vector field X € C°^{M, TM) 



(3.4) ^(pi,P2:P3,P4) ~ 

P2 

163 In order for this transformation to be invariant under scaling, the transformation is nor- 

164 malized in [22] by multiplying each face normal with the inverse of the square root of its 

165 norm. The normalized vector field factors through the quotient M / ~. We propose a differ- 

166 ent normalization, which also preserves the direction of vectors on M and the singularities 

167 of the vector field, which makes dynamical studies simpler. 



1.(4,3,2) 
1.(4,1,3) 
^.(4,2,1) 
^.(1,2,3)1 



for the tetrahedron 
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168 Definition 3.2. For a tangent vector v on M let 

169 The scaling invariance is formulated more formally as follows. 

170 Lemma 3.3. Let X be a homogeneous quadratic vector field on M , that is, the coefficients in 

171 the canonical basis are homogeneous quadratic polynomials with respect to pi . Furthermore, 

172 let X be a translation invariant. Then Y = '^{X) is a vector field on the quotient M/ ^. 

173 More precisely, the vector field Y on N is given by Yq = I?(7r)(^(Xp)), where Tr{p) — q and 

174 Z5(7r) : TpM T^S'^^""!)"! is the differential ofn: Af -> 53(n-i)-i^ 

Proof. Let X' = 'J'(X) and notice, that tX' = XL and X' = X' , y We need 



to check that -D(7r)(Xp) = D{iT){X'p,) for p ^ p' . Clearly, if p' = p + {po, ■ ■ ■ ,Po), then 
D(7r)(X;,) = -D(7r)(X;). If p' = Xp, then 

D{7r){X'^,) = D(a){X[{\p) - X[^{\p), . . .,X'^_,{Xp) - Xl^Xp)) 

= j^a{Xp + t{X[{Xp) - X^Xp), X:,_,{Xp) - X'„{Xpmt=o 

= j^a{X{p + t{X[{p) - XUp), ■ . ■,X:,_,{p) - X'M)))\t=o 

= ^a{p + tiX[{p) - X'M, ■ ■ -^K-M - ^;(p)))Uo = D{t:){XL). 



175 □ 

176 The remaining work rests on a generalization of the following simple observation to 

177 polyhedra and polyhedral meshes. If p is a tetrahedron mesh, then the gradient field of its 

178 oriented volume 

(3.5) VOl(p) = i((p2 -Pl) X (P3 -Pi)) • {Pi-Pl) 

179 is 

ix = Vvol, 
6 

180 where X is given by Equation (3.4). The dynamical behavior of the geometric element trans- 

181 formation for tetrahedra seen in [22] therefore has a classical answer; it is the discretization 

182 of the gradient flow of the volume function on M . In order to generalize this, we need to 

183 discuss gradient fields in a more general context. 

184 4 Gradient fields 

185 The gradient flow is an important concept and a powerful tool with many applications in 

186 topology, global analysis, dynamical systems and mathematical physics. Perelman observed 

187 in his proof of Thurston's geometrization conjecture [10], that one can interpret Hamilton's 

188 Ricci flow as a gradient flow. Morse theory [14, 3] with all its generalizations and infinite- 

189 dimensional manifestations in mathematical physics like Chern-Simons theory [23], Yang- 

190 Mills theory [1, 2, 15] and other gauge theories deduces topological information from the 

191 study of the singularities and the flow of a gradient field. It is a natural consequence of the 

192 strong results surrounding gradient fields and their flows, that we should make use of them 

193 to shed some light on the dynamics of GETMe and generalize it. In this section we will 

194 discuss some important properties of gradient fields. In particular we will show, that the 

195 (normalized) gradient flow starts and ends at singularities of the gradient fleld, and we will 

196 characterize these singularities by a simple equation. 

197 We denote \yy X ■ p the real- valued function f^{p) ■= Xp ■ p on M. Notice that vol is of 

198 the form f-^. Even though / is not a function on the quotient A// ^, / restricts to iV C M 

199 and we can consider its flow curves (or integral curves), each of which depends on a time 

200 parameter t € R. Since A*" is a compact manifold, all vector fields on N are complete, i.e. 

201 their flow curves exist for all time. 
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Proposition 4.1. Let X he a translation invariant vector field on M . Then the flow of the 
vector field Y =^ D{tt){X) on N d M has a singularity at p d N if and only if 

(4.1) Xp = T{Xp) for some X — X{p). 

Proof. For p Cz N d M we have cr(p) = t{p) — p. A singularity p (z N of Y is characterized 

by 

O^D{n){Xp)^D{a)oD{T){Xp). 

The translation invariance of X and the identification of the tangent space at any p with 
R^" implies that D{T){Xp) — T{Xp), and we can rewrite the above as 

= D{a)r(^p^riXp). 

Let us identify Im(L)(T)p) =R3("-i). Then the kernel of i:>((T)p|c(^)^ : R3("-i) TpS^'^'^-^^'^ 
consists of the radial vectors, i.e. vectors of the form Xp. Therefore, Yp = if and only if 
T{Xp) — Xp for some A = X{p). □ 

The following fundamental quality of gradient fields is the main ingredient in the regu- 
larizing behavior of the X or its normalization. 

Theorem 4.2. Consider the function f{p) := X ■ p: A/ — > R on (Af, (•, •)), such that the 
components of the vector field V/ are translation invariant quadratic homogeneous vector 
fields. Then Y = D{Tr)(\/f) is the gradient field of /Iat, and the flow lines of Y and 
Y = _D(7r)(5'V/) on N <Z M start and end at singularities of Y . 

Proof. Let 7 be a flow line of Y on iV C M . Since V/p ~ frip) we have for p E N 

D{T)Vfp = Vfp. 

Since D{a) is given by orthogonal projection on T^^^-^N and dividing by we have 
4.2 D{TT)Vfp - . 

In particularly we have D{'k)V fp = V fp — fp ■ p) ■ p for p E N. Therefore, for Zp e TpN 

(V(/U)p,Zp)^ = D{f\N){Zp) = DfiZp) = {Vfp,Zp) = {D{^)[Vfp) + {Vfp-p) -p.Zp) 
= {Di7r){Vfp),Zp) = {D{7r)iVfp),Zp)h. 

This shows that D{tt){V fp) is the gradient field of /|jv- Then 

|(/o7) = (i^(7r)(V/,(,)),7(t)>. = ||7(<)ll^. > 

implies that / o 7 is increasing for all t. 

Since N is compact and / continuous, Im(/ o 7) is a bounded set, and thus 

.^"^ :^(/°7) = o. 

t-i.±oo dt 

Since N is compact, we can find a sequence i„ of real numbers such that 

lim t„ — 00 and lim 7(t„) ~ p E N. 

n— fC30 n— fcjo 

By the continuity of / we get Yp = and p is singular. A similar argument shows that 7 
starts at a singularity. 

Similarly to Equation (4.2) we have L)(7r)(\I'V/p) f-„{p) for p E N. For a flow line 
7 of y we therefore have 

|(/°7) = (V/^(t),7(i)) = (V/^(,),D(^)(vE'V/^(,))) = (V/^(,),vl>V/,(,)) > 0. 
The rest of the argument remains the same. □ 
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228 We are interested in the normalization of X, because '^{X) is a scaling invariant 

229 vector field if X is quadratic homogeneous. One can also ask about the dynamic behavior of 

230 normalizations given by division of other powers of the norm and coordinate- wise normaliza- 

231 tions. We expect them to have the same qualitative behavior. In the case of simply dividing 

232 by the norm of X , this dynamical behavior is described in Thom's gradient conjecture [12]. 

233 Theorem 4.2 suggests that X and "i^iX) are topologically equivalent, that is, there is a 

234 homeomorphism N N carrying trajectories to trajectories and preserving the direction 

235 of increasing time. In addition, the singularities of X and "^{X) are the same. For our 

236 purpose it is sufficient, that X and "^(X) have the same qualitative behavior as described 

237 Theorem 4.2. Therefore, we only need to analyze the flow of X. For a function /: Af — > R 

238 we let Vj/ :— {dp.^f.dp.^f.dp.^f) be the gradient of / with respect to pj only. We leave 

239 the straight-forward proof of the following result to the reader. 

Lemma 4.3. For p ~ {pi, . . . ,p„ ) e M we have 

Vi((pi X P2) - Pa) =P2 X P3, 
V2((pi X P2) -ps) =P3 X pi, 
V3((Pl X P2) -pz) =Pl X P2, 

{Pl X P2) ■ P2 = (P2 X Pi) • P2 = 0, 

71 

and (V[(pi X P2) ■ p^]) ■ P {^i[{Pi x P2) ■ Ps]) ■ Pi = 3(pi x P2) ■ Ps- 

i=l 

240 The above lemma allows us to characterize vector fields Y such that V/"^ = Y . 

241 Proposition 4.4. Let X he a vector field on M , such that Xj is a linear combination of 

242 cross products in the variables pi, , j,i — 1, . . . , n. Then for Y ~ V/"^ we have 

v/^ = sr. 

243 Proof. It follows immediately from Lemma 4.3 that V/"^ • p = 3X ■ p. Therefore, 

v/^ = v(v/^ . p) = v(3(x • p)) = 3V(Xp • p) = sr. 

244 □ 

245 We have already seen that the vector field X for tetrahedra given in Equation (3.4) is a 

246 gradient field on M. In order to generalize the construction of gradient fields to arbitrary 

247 polyhedra and meshes, we need to discuss triangulations. 



5 Triangulations 



249 The notion of a triangulation exists in all dimensions. A triangulation of a (convex) poly- 

250 hedron is a set of tetrahedra, such that the union of all of them equals the polyhedron, 

251 and the intersection of any pair of them is a (possibly empty) common face. It is most 

252 convenient for our purpose to consider only triangulations such that the union of the ver- 

253 tices of all tetrahedra agrees with the vertices of the polyhedron. This can be generalized 

254 to polyhedron meshes, where we require the triangulation of the mesh to be induced by 

255 a triangulation of the underlying polyhedron. More formally, a triangulation of a polyhe- 

256 dron mesh (C,p) is a set of tetrahedron submeshes {(C(P'),p*)}i=i_...,fc with IJj = |C|, 

257 V)i^{C{P^)) = y{C) and p' = p\v{c(P^))- element-wise triangulation of a polyhedral 

258 mesh {[jiC{P^),p) is given by T = Ui where Ti is a triangulation of (C(P*),p'). If we 

259 fix an isomorphism class of a polyhedral mesh and an identification with (R'^)", then an 

260 element- wise triangulation of p = (pi, . . . ,p„) can be written as a set of tetrahedron meshes 

261 p* of p, which are 4-subtuples of p. For example a pyramid mesh p = (pi,P2,P3,P4,P5) as in 

262 Figure 2 can be triangulated in two different ways Ti = {(pi,p2,P3,P5), (pi,P3,P4,P5)} and 

263 T2 = {(P1,P2,P4,P5), (P2,P3,P4,P5)}- 

264 This allows us to construct a transformation or equivalently a vector field for an arbi- 

265 trary (polyhedral) mesh from the vector fields associated to the tetrahedra of element-wise 

266 triangulations. 
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Figure 2: The two different triangulations of a pyramid 



Definition 5.1. Fix the isomorphism class of a polyhedral mesh and an identification with 
(R'^)". Let T be a subset of all element- wise triangulations of p and |T| the cardinality of 
T. Then we let 



J2Xp. forT = {p\...,/-}er 



and define the field 



X 



T 



where is given by X in Equation (3.4) applied to T. 

This vector field will yield a flow on the manifold of polyhedra and polyhedral meshes. 
We will see in Sections 6 and 7 two which extent the regularizing behavior for tetrahedra 
can be observed and proven for more general polyhedra. 

Let us compute for the polyhedron whose triangulations consist only of two tetra- 
hedra: the pyramid. The two triangulations in the set of all triangulations T are shown 
in Figure 2. If X^, i = 1, . . . , 4, p C {pi, . . . ,^5} are the gradient fields associated to the 
four tetrahedra shown in the two triangulations, the vector field for the pyramid is given 
by = lEti^'- By linearity 3X^ = ^ Eti = V(ix; -p) - ViXj ■ p), and 
therefore X^ is a gradient field. We compute 



(5.1) 



xT = 



1^(5, 4,2)+ K5, 4,3, 2)- 
:^(5,1,3) -hJ/(5,l,4,3) 
i.(5,2,4) + K5,2,l,4) 
1^(5,3,1) + i^(5, 3, 2,1) 
2.1.(1,2,3,4) 



It is also straight-forward to construct the vector field for the prism, because the set T 
of all triangulations has only six elements. 



(5.2) 



xj 



1/(3, 2, 4) +1.(2,5,4, 6, 3)' 
1/(1,3,5) +1.(3,6,5,4,1) 
i/(2,l,6)+Kl,4,6,5,2) 
:/(5,6,l) + (.(6,3,1,2,5) 
K6,4,2) + K4,l,2,3,6) 
1/(4,5,3) + i.(5, 2, 3, 1,4) 



for the prism 




We have observed that the vector field X for tetrahedra as in Equation (3.4) is the 
gradient of the volume function (3.5). Clearly, for triangulations of convex polyhedra the 
volume function is additive: 

k 

vol(P)=^vol(P'), 

where T — {P*} is a triangulation of the polyhedron P. Therefore we define for an element- 
wise triangulation T = IJ^ e T of a polyhedral mesh p 



voF(p) := 
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288 We illustrate the dependence on the element-wise triangulation T by way of the pyramid 

289 mesh in Figure 3, where we omit some faces in order to clarify the location in three dimen- 

290 sional Euclidean space. Figure 3 shows a constellation, for which the two triangulations 
from Figure 2 give non-zero volumes with opposite sign. 




Figure 3: A pyramid mesh 

291 

292 Due to this dependence we can define a useful volume function by averaging over all 

293 element-wise triangulations T G T. 

294 Definition 5.2. Let p be a polyhedral mesh, then we define its mean volume to be 
(5.3) vol:--^^voF. 

295 The mean volume for meshes generalizes the (generalized) volume for triangulated poly- 

296 hedral surfaces introduced by Connelly [6], which was used in the proof of the bellows 

297 conjecture [5, 16]. Note that we can equivalently express the mean volume as a sum of 

298 <lei{pi,pj,pk) over all triangles (j>i,Pj,Pk) in all triangulations of the boundary surface. The 

299 contribution of quadrilateral (non-planar) facets to the mean volume can be reformulated 

300 using certain doubly-ruled surfaces [7]. In applications, we usually only encounter meshes 

301 with triangular and quadrilateral facets, and the use of the doubly-ruled surfaces allows us 

302 to treat the polyhedron meshes as 3-dimensional objects. Let us summarize the preceding 

303 elementary discussion in the following result. 

304 Theorem 5.3. LetT = [J^Ti be an element-wise triangulation of a polyhedral mesh {(C{P^),p)}i^i^,,,^k, 

305 where Ti is a triangulation of {C{P^),p^). Then the gradient of 

k k 

voF(p) = ^ voF(p*) = E E 

4=1 i=l qeTi 

306 is given by 

VvoF = -X'^. 
6 

307 Remark 5.4. It is not difficult to see, that the gradient of the mean volume Vvol preserves 

308 the centroid 1/n ■ "^pi of a given polyhedron p — (pi, . . . ,Pn)- This is an important fact 

309 for mesh smoothing applications, but in our work this does not play a role, so we will not 

310 discuss this any further. 

311 We can construct concrete gradient vector fields for other polyhedra more easily from 

312 symmetry-invariant subsets of triangulations. This can be seen for the hexahedron, which 

313 has a two-element symmetry-invariant subset. Let us therefore formally introduce and dis- 

314 cuss symmetry. 

315 6 Symmetry 

316 Mani [13] showed that for each polyhedral graph G, there exists a convex polyhedron Pq 

317 such that every automorphism of G is induced by an isometry of Euclidean space. We call 

318 Pg a symmetric polyhedron. Examples for symmetric polyhedra include the platonic solids, 

319 semi-regular polyhedra, the right regular prism, right regular anti-prism and the regular 

320 pyramid. In this section every polyhedron P is assumed to be symmetric. 
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Consider the symmetry group F of a (symmetric) polyhedron P C as symmetries of 
R^, that fix P. Therefore g £ T acts on all points p E R^ as usual hy g-p G R'^. In particular, 
it naturally acts on a triangulation {P*} of P, and we can consider a F-invariant subset of 
triangulations of polyhedra. Similarly, it also acts on a triangulation {(C(P*),p')}i=i_...^fc of 
a polyhedron mesh via 



5-{(C(n,pO},.i,...,fc := mg ■ P^),f)h 



=i,....fc- 



In this way, a F-invariant triangulation of a polyhedron gives a F~invariant triangulation 
of a polyhedron mesh. Notice, that the orbit of an arbitrary triangulation under F yields a 
F-invariant subset of all triangulations. 

Let us consider the two possible triangulations of the hexahedron consisting of 5 tetra- 
hedra each and apply this theorem to this F-invariant subset T of all triangulations. 



(6.1) 



xj 



1/(2, 5, 4) + 1/(6, 5, 8, 4, 3, 2)' 
i/(3,6,l) + z/(7,6,5,l,4,3) 
1/(4, 7, 2) + 7, 6, 2, 1,4) 
//(1,8,3) + i/(5,8,7,3,2,l) 
Kl, 6, 8) + 1/(6,7,8,4, 1,2) 
1/(2,7, 5) + 1/(7, 8, 5, 1,2, 3) 
1/(3, 8, 6) + 1/(8, 5, 6, 2, 3, 4) 
1/(4, 5, 7) + 1/(5, 6, 7, 3, 4,1) 



P8^ 



P5 



for 



P6 



P7 



^P3 



Pi P2 



In addition to the hexahedron, it is now a straight-forward task to write down and check 
gradient fields for the rest of the platonic solids. The gradient field for the octahedron can be 
found in Section 7.3. In Section 7.6 we mention a gradient field for the icosahedron, which 
is a variation of the gradient field of the mean volume. We leave the explicit gradient field 
for the dodecahedron to the interested reader, because the expression is considerably longer. 
The symmetry group of P also acts on the vertices Vi , . . . , of P by permuting them, that 
is, Vo-g(i) '■= g ■ Vi for some ag G 5„. Then F acts on the mesh p e (R^)" of C{P) from the 
left as usual g-p= (g ■ pi, . . . , g ■ Pn) as weh as from the right viap-^ = (Pa^ii), ■ ■ ■ ,P<Tg(«))- 
We call a mesh p symmetric, ii p — g ■ p ■ g^^. 



7 Regularity 

Intuitively, symmetry and regularity are related notions. One definition found in the litera- 
ture is, that a polyhedron is regular, if the action of its symmetry group acts transitively on 
its flags, so that the platonic solids are the only regular, convex polyhedra [8]. Nevertheless, 
prisms and pyramids are also sometimes called regular, if they are symmetric. This notion 
can also be extended to isogonal polytopes. Let us therefore call a pyramid, a prism and an 
isogonal polyhedron regular, if it is symmetric (and convex). It has been observed heuristi- 
cally in [19] , that the iterative application of certain actions induced by the symmetry group 
on polygons regularizes them. In a way this is our hope for the flow of the vector fields X'^ 
for a F-invariant set of triangulations T- In this section, we will discuss this in more detail. 

In Definition 5.1 we have defined a vector field, which is the gradient field of the mean 
volume function. In the case of tetrahedra, we will show directly in Section 7.1 that the 
non-coUinear singularities are regular tetrahedra. In particular we see by Table 1, that the 
positively oriented one is a sink, and the negatively oriented one is a source of the gradient 
flow. 

It is tempting to think that this generalizes to other polyhedra with non- vanishing mean 
volume, namely that the singularities of the mean volume flow on N are regular polyhedra in 
some sense. Unfortunately, this naive approach will only work for rather simple polyhedra, 
even if they have a large symmetry group. For example, in the case of the icosahedron, there 
are a lot more singularities of the mean volume function than just the regular icosahedron, 
and, in particular, the regular icosahedron does not maximize the mean volume function 
on N. In this case, a slightly more complicated gradient fleld seems to have the desired 
symmetric icosahedra as singularities: the regular icosahedron and the great icosahedron 
(one of the Kepler-Poinsot polyhedra). 
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364 While we have confirmed numerically, that the gradient flow of the mean volume exhibits 

365 a regularizing dynamic behavior for the tetrahedron, the pyramid, the octahedron, the 

366 hexahedron and even the dodecahedron, we can rigorously confirm this behavior only for 

367 the tetrahedron, the pyramid and the octahedron. For the hexahedron, we can show it for 

368 a slightly different gradient field. Surprisingly and regrettably, we have not yet been able 

369 to prove it for the prism, even though numerical tests suggest, in particular a computation 

370 of the eigenvalues of the Hessian of the coUinear prisms, that the mean volume is a Morse- 

371 Bott function just like we prove it for tetrahedra in Section 8. The icosahedron and the 

372 dodecahedron have only been studied numerically. 

373 Equation (4.1) for A 7^ can be reformulated more geometrically. For any singularity, 

374 there exists pq G R"^ satisfying Xp — (po, • • ■ ,Po) + -^P- For A ^ and 

1 = j{Po, ■ ■ ■ ,Pq)+P 

375 this is equivalent to Xq — Xp — Xq. In other words, for Xp ^ {po, . . . ,po) singular p are 

376 being rescaled via Xp. We make the following definition: 

377 Definition 7.1. Fix an isomorphism class of a polyhedral mesh {C,p) and let X be a 

378 gradient vector field on N — M / ~ given by a linear combination of cross products. The 

379 mesh p is X -optimal, if there exists q ^ p satisfying 

(7.1) Xq = Xq for some A 7^ 0. 

380 In particular, no coUinear p is X-optimal. The sign of A determines the orientation of p. 

381 We allow Xp' = for polyhedron submeshes p' of q as long as Xq 7^ 0. 

382 We expect that the X-optimal meshes are sinks or sources depending on the sign of A. 

383 It is much more difficult to determine the shape of the singularities. Let us therefore discuss 

384 X-optimality for the platonic solids, the prism and the pyramid, when X is the gradient of 

385 the mean volume or a slight variation of it. 

386 7.1 The tetrahedron 

387 The tetrahedron is the simplest case. We show in Section 8, that the mean volume is a 

388 Morse-Bott function on the sphere N. In particular, the dynamics of its gradient flow are 

389 well-understood, see Section 9. 

390 Theorem 7.2. The non-collinear singularities of the gradient of the (mean) volume for 

391 tetrahedra are regular tetrahedra. 

392 Proof. Let X be the gradient of the volume given in (3.4). If p is X-optimal, let us assume 

393 without loss of generality that p4 = 0. Then 

^{Pi -P2) = ^(^(4,3,2) - i.(4,l,3)) =P3 X Q(pi +P2)) . 

394 Therefore pi —p2 is orthogonal to the hyperplane spanned by ps and 1 {pi +p2 ) . In particular, 

395 the face with vertices pi, p2 and ps is isosceles, see Figure 4. By symmetry, all edges of the 
tetrahedron are of equal length, and therefore the tetrahedron is regular. □ 



P3 




Pi P2 



Figure 4: Isosceles triangle 

396 
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7.2 The pyramid 



398 
399 
400 
401 



Let X be the gradient of the mean volume function for pyramids given in (5.1). Suppose p 
is X-optimal and assume P5 = without loss of generality. We compute P4 — = P2 — Pi 
for every X-optimal p £ N. Therefore the base of a singularity is a (planar) parallelogram. 
Furthermore we find for A 7^ as in Definition 7.1 

^(P2+P4) = 3^.(4,3,2,1) and ^(pi +^3) = 3^.(4, 3, 2, 1). 
Therefore ps lies vertically above the center of the parallelogram base. Now consider 

A(P4 - Pi) = 1^(5, 3, 1) + ^.(5, 3, 2, 1) - i.(5, 4, 2) - i/(5, 4, 3, 2) 
= 1/(5, 3, 1) + i/(5, 2, 4) + j/(5, 2, 1) + i/(5, 3, 4). 

We observe that for some constant C 



iy{5,3,l) + 1^(5,2,4) = C{p4- Pi) 



404 
405 



and 

(K5,2,l) + K5,3,4))±p2-Pi. 

Since p is non-collinear, the base is a parallelogram and does not lie in the base, we have 
p{5, 2, 1) + ^(5, 3, 4) ^ 0. We conclude that p4 —pi and p2 ~Pi are non-zero and orthogonal 
to each other. We have using p^ — 

A(P4 - P2) = 2^/(5, 3, 1) + i^{5, 3, 2, 1) + ^.(5, 3, 4, 1) 

= 2p3 Xpi + (p3 - pi) X {p2 + P4) 
= (P3 - Pi) X (pi + P3 +P2+ Pi). 

The right-hand side is orthogonal to p^ — pi- Therefore, the base is a square, and p is a 
symmetric pyramid, whose height is determined by the following result. 

Theorem 7.3. Every X -optimal pyramid is of the form 7r(p) G N , where 

(0,0,0) ■ 

(2,0,0) 
(2,2,0) 
(0,2,0) 
(1,1,75). 

Proof. Consider a pyramid of height h with square base of area and apex orthogonally 
above the middle of the base. If p is X-optimal, then we compute 2, 3, 4)|| = 2a^ and 
||z/(5, 2, 4)11 — h ■ \\p4 —p2\\. For A 7^ as in Definition 7.1 we have 



2Xh = 2A 



P1+P2+ P3 + Pi 
4 



= 112^.(1,2,3,4) -^^4, 3, 2,1)11 =5a2 



and 2A||p3-pi|| = ||4i/(5, 2, 4) + ^2, 1, 4) + i.(3, 4, 2)|| - 4/i||p4 - P2II 



Since ||p3 ~ Pi\\ = \\pi — P2\\, we compute 



h V5 



□ 



408 
409 
410 
411 
412 
413 
414 
415 
416 



Remark 7.4. From the perspective of the finite element method, it is desirable, that all 
edges have the same length. We have two options to take care of this issue. We can ignore 
it entirely by arguing, that the edges connected with the apex for p as in Theorem 7.3 have 
length \/7 « 2.64575, which might be close enough to 2 for mesh smoothing applications. 
Alternatively, if we are willing to let go of the theoretical advantages of a gradient field, it 
is straight-forward to vary the vectors in the gradient field so that all edges have the same 
length and the optimal singularities are preserved. Clearly, we have a lot of (topologically) 
equivalent transformations at our disposal and we can control the shape of the optimal 
pyramid, if wc could only prove topological equivalence of the corresponding vector fields. 
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417 



7.3 The octahedron 



419 
420 
421 
422 
423 
424 



425 
426 



It is straight-forward to compute the gradient field of the mean volume to be 

Pi/ 



(7.2) 



1 



K2,3,4,5) 
Kl,5,6,3) 
Kl,2,6,4) 
Kl,3,6,5) 
Kl,4,6,2) 
K2,5,4,3) 



for P2 




Let p be X-optimal. Then we can readily calculate that Ps — Pe = Pi ~ P5- Together with 
analogous statements for the other edges, we sec that C(2,3,4, 5) and the other equator 
curves are a parallelograms. Furthermore, pi — pe = 2i/(2, 3, 4, 5), therefore pi ~ pe is 
orthogonal to the plane, in which C(2,3,4, 5) lies. Together with the analogous statements 
for the other cases this implies orthogonality of the planes, in which C(2, 3, 4, 5), C(l, 2, 6, 4) 
and C(l,3,6,5) lie. 

To prove the regularity of the octahedron, it remains to check that the area of the rhombi 
circumscribed by C(2,3,4, 5), C(l,2,6,4) and C(l,3,6,5) are all equal. We compute 

11^/(1,3,6,5)11 = -Pell • \\P3-P5\\ = ^1121.(2,3,4,5)11 • 11^1,2,6,4) - ^^(1, 4, 6, 2)|1 
= 2A-2|1j.(2,3,4,5)|1- 11^/(1,2,6,4)11. 

By symmetry we also have ||:/(2, 3, 4, 5)|| = 2A-2||;y(i, 3, 6, 5)|| • ||!/(1, 2, 6, 4)|| . Therefore, 
once again by symmetry, 

11^/(2,3,4,5)11 = |li/(l,3,6,5)|l = |li/(l, 2, 6, 4)|1 - 
This shows the following. 

Theorem 7.5. Every X~optimal p is a regular octahedron. 



429 7.4 The prism 

430 We are in a bit of a dilemma in the case of the prism. We have not been able to determine 

431 the singularities of the gradient X of the mean volume on N given in Equation (5.2). The 

432 computation of the optimal singularities for the simpler vector field 



z/(3,2,5,4,6y 
i/(l,3,6,5,4) 
1/(2,1,4,6,5) 
1.(5,6,3,1,2) 
1/(6,4,1,2,3) 
1/(4,5,2,3,1) 



433 is straight-forward. Let p be K-optimal. Then we have Ps ^ Pe = P2 ^ P5 and analogous 

434 equations for the other quadrilateral facets. This proves that the quadrilateral facets are 

435 parallelograms. Therefore i/(l,2, 3) = — 1/(4, 5,6) and pi — p4 = i/(l,2,3) — z/(4, 5,6) = 

436 2i/(l,2,3). Every K-optimal p is therefore a symmetric prism of a certain height. Our 

437 dilemma is that we can only prove half of the facts necessary for a complete picture of the 

438 dynamics of X and Y: Unfortunately Y is not a gradient vector field, as can be checked 

439 by Proposition 4.4, but X is; on the other hand, we know the singularities of Y , but not of 

440 X. Based on numerical tests we conjecture, that Y is topologically equivalent to X. This 

441 would solve our dilemma. 

442 Assuming that X and Y have the same number of optimal singularities, the following 

443 result determines all X-optimal p. 



(7.3) 



Yr, 
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Theorem 7.6. Every symmetric X~optimal prism is of the form t:{p) G N , where 



445 
446 
447 
448 



450 
451 
452 
453 



(0,0,0) 
(2,0,0) 
(1,V3,0) 

(0,0, 
(2,0, 
(1,73, 



Proof. Consider a prism with equilateral triangles for the base faces and rectangles for the 
other three faces, where p G N. Let a be the length of each side of the bases, and let h be 
the height of the prism. Therefore the area of each base is \/3a^/4, and the area of each 
rectangular face is h ■ a. The prism p is optimal, if Xp = t(X). Therefore 



(7.4) 



and 



/ TT 

\a = X\\p3 -P2II = ||^p3 - XpJ\ = 2i7cos l^- 



hV3,., 



where H = -h ■ a 
2 



\h = X\\p4 ~pi\\ = \\Xp^ - Xp^W = a^\/3. 



The argument behind Equation (7.4) is indicated in Figure 5; an elementary calculation 
shows that the length 2H of the projections of 2Xp^ = v{2, 1, 6) + z^(l, 4, 6, 5, 2) and 2X^2 = 
v{l, 3, 5) + 1/(3, 6, 5, 4, 1) to the base each agrees with 3 times the area of each rectangular 
face. It follows that 

h fo. 

□ 




Figure 5: The projection of the prism to the base 

454 



455 Remark 7.7. As in the case of the pyramid, a prism with all edges of the same length might 

456 be more desirable. It is straight-forward to construct a non-gradient field from the gradient 

457 field, such that the optimal singularities have the desired edge-length. 

458 7.5 The hexahedron 

459 We have not been able to compute the optimal polyhedra for the gradient of the mean 

460 volume given in Equation (6.1). However, the gradient vector field 



1/(3, 6, 8) + v{2, 5, 4) + j/(6, 5, 8, 4, 3, 2) 
1/(4, 7, 5) + z/(3, 6, 1) -t- 1/(7, 6, 5, 1, 4, 3) 
1/(1, 8, 6) + j/(4, 7, 2) + 1/(8, 7, 6, 2, 1, 4) 
1/(2, 5, 7) + v{l, 8, 3) + 1/(5, 8, 7, 3, 2, 1) 
1/(2, 7, 4) + 6, 8) + 1/(6, 7, 8, 4, 1, 2) 
1/(3, 8, 1) + v{2, 7, 5) -I- 1/(7, 8, 5, 1, 2, 3) 
1/(4, 5, 2) + z/(3, 8, 6) + 1/(8, 5, 6, 2, 3, 4) 
v{l, 6, 3) + v{4., 5, 7) + t/(5, 6, 7, 3, 4, 1) 



(7.5) 



^ 2 



14 



February 26, 2013 



461 differs from the gradient of the mean volume only slightly. In order to compute the optimal 

462 singularities, let us denote the faces by Ai and the corresponding face normals hy rii, i = 
1, . . . , 8 given by the cross product of the diagonals as indicated in Figure 6. 

.n4 




Figure 6: The face normals of the hexahedron 

463 

464 It is straight-forward to compute pi — pi = P5 — Ps for any y-optimal polyhedron p. 

465 It follows that pi — Ps = (pi — Ps) + {P5 — Ps), so that A3 is planar. By symmetry every 

466 singularity has planar faces. It also follows frompi— p4 = P5—ps that yli|j>li_|_3 for i = 1, 2, 3. 

467 Furthermore, we have A(p4—p8) = fti +i^(2, 5, 7) — ri4 — 1^(1, 6, 3). Notice that (i^(2,5,7) — 

468 1/(1, 6 , 3))||ni, so that {p4—ps)\\ni and the faces A3 and Ai are orthogonal to each other. By 

469 symmetry the same is true for the other faces. Such an element is known as an orthotope. 

470 For an orthotope it is straight-forward to see that 

^(ni +n2 +713) = 1/(3,6,8). 

471 Therefore, the first entry of Yp satisfies 

1/(6,5,8,4,3,2) = ^(ni+n2 + n3) + z/(2,5,4) = K2,5,4) + i/(3,6,8), 

472 and the analogous relation holds for the other entries. Let us compute the area of ^3 of our 

473 orthotope in two equivalent ways. On the one hand the aria is simply ||n3||, on the other 

474 hand it is \\pi — Ps\\ • \\P5 ^PsW- It follows that 

ll^sll = 4A^||ni - ^411 • \\n2 - nsll = 16A^||ni|| • ||n2||. 

475 By symmetry we have ||ni|| = 16A^||n2|| • ||n3|| and therefore ||n3|| — 16^A'*||n2||^||ri3||. We 

476 compute ||n2|| = (16A)^^. By symmetry we see that ||rii|| = ||7^-2|| = ll^sll — (16A)^^. 

477 Therefore every optimal singularity is a cube. 

478 7.6 The icosahedron and the dodecahedron 

479 Due to its triangular faces the i-th entry of the gradient X of the mean volume is given by 

480 the normal vector to the surfaces enclosed by the link of pi. However, the singularities of 

481 the gradient fiow are not only regular icosahedra. In fact, the regular icosahedra are not 

482 maximal with respect to the volume function. We can vary the gradient fiow by adding 

483 the normal vector to the surfaces enclosed by the link of the point opposite to Pi. More 

484 precisely, the resulting gradient vector Y is given by Equation (7.6), where only the first 

485 entry is shown and the other entries can be found by rotation symmetry. Experiments have 

486 shown that the y-optimal singularities are regular icosahedra and regular great icosahedra. 



(7.6) 



X 



r 



K2,...,6) 
+K7,...,11) 

-j/(5, 7, 3, 10, 6, 8, 4, 11,2,9) 



for 




P12 
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487 If X is the gradient of the mean volume for the dodecahedra, we have observed experi- 

488 mentally, that the X-optimal singularities are regular dodecahedra. 

489 8 Volume as a Morse-Bott function 

490 In this section we show that the signed volume is a Morse-Bott function on the 8-sphere of 

491 tetrahedron meshes. In particular this explains the qualitative behavior of its flow. 

492 In general, the group S0(3) acts on (pi, . . . ,p„) g M by rotating all pi simultaneously, 

493 and the vector field V vol is invariant under rotation. Therefore, all critical points of the 

494 mean volume on N are degenerate due to the rotational invariance. The action of SO (3) on 

495 the (3n — 4)-sphere N is proper and free as long as the pi are not collinear. The stabilizer 

496 Sp = {A G S0(3) \p - A — p} ioi every collinear p G M is . The image of the submanifold 

497 of collinear p G M under r: M — > R"~^ is diffeomorphic to x R"^^\{0}, where the 

498 first factor stems from the rotation of the collinear rays identified with the second factor. 

499 Therefore, the subspace Sc of collinear p G TV is diffeomorphic to x S'"~^ via a o t. 

500 Clearly, the set of singular tetrahedra with vanishing volume consists of the collinear ones 

501 and is diffeomorphic to x S^. In general, Sc is only a subset of the critical points of the 

502 mean volume. 



In order to show that the (mean) volume on N is Morse-Bott for tetrahedra, we need to 
compute the eigenvalues of the Hessian. The Hessian of the volume on M is equal to the 
gradient of V vol. In order to compute the eigenvalues of the Hessian of the restriction of 
/ : M -> R to at a regular polyhedron p, we can compute the eigenvalues of the Hessian 
of /oTT and keep in mind, that six eigenvalues are zero; three from the translation invariance, 
two from the rotation invariance and one from the scaling invariance. On the level of vector 
fields we need to take the push-forward TT^,{Xp) of Xp via tt* = (cr o r)* : TpM — >■ TTr(p)N 
before we take the gradient. As we have seen earlier, 

T* {Xi, . . . , Xn) = {Xi — Xn, . • . , Xn~l ~ Xn, 0) 

and a* {Xp) — p — '^^^ ■ Xp, 



503 where we assume that N = N x {0} C M. In order to compute the eigenvalues of the 

504 Hessian of f\j\j at positively oriented polyhedra, we can compute the eigenvalues of the 

505 Jacobian matrix given by the gradient of TT^,{Xp) C TpM. The result of this computation is 

506 summarized in Table 1. By symmetry, the eigenvalues for the negatively oriented polyhedra 

507 are the same, only with the opposite sign. This agrees with the observation that the mean 

508 volume has a global minimum at the negatively oriented regular polyhedron and a global 

509 maximum at the positively oriented regular polyhedron. 



Type 


X 


n 


(A,m) = eigenvalue A of tt^{\/X) with multiplicity m 


tetrahedron 


(3.4) 


4 


3 


(-^8/3, 6) 


pyramid 


(5.1) 


5 


3 


(-V20/7,6),(-y577,3) 


octahedron 


(7.2) 


6 


3 


(-4/^3,6), (-2/^3, 6) 


prism 


(5.2) 


6 


3 


(-V3, 6), (-2/V3, 2), i-V^/2, 2), (-1/^3, 2) 


hexahedron 


(6.1) 


8 


3 


(-73, 6), (-5/712, 1), (-2/73, 3), (-73/2, 3), (-1/73, 5) 


hexahedron 


(7.5) 


8 


3 


(-4/73, 6), (-2/73, 12) 



Table 1: The eigenvalues of the Hessian at 
positively oriented regular polyhedra 



510 We also compute that the Hessian at collinear tetrahedra has two positive and two neg- 

511 ative eigenvalues. This shows that the (mean) volume is a Morse-Bott function on the 

512 8-sphere N, more precisely, the critical sets are submanifolds S^, S"^ x S"^ and S"^ at the 

513 maximum, level and the minimum, respectively, and the Hessian is non-degenerate in 

514 the normal direction. Furthermore, the positively oriented regular tetrahedron is a sink of 

515 the gradient flow as well as of the flow of \l/(Vvol), while the stable and unstable mani- 

516 fold are two-dimensional for collinear tetrahedra. Generically, all tetrahedra are therefore 

517 regularized by the volume gradient flow, more precisely, the set of all tetrahedra, which 
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518 are regularized, is open and dense in N. Alternatively, a small perturbation of the mean 

519 volume in a neighborhood of the level singularities is a Morse function with two singulari- 

520 ties consisting of the positively and negatively oriented tetrahedra. Unfortunately, for other 

521 types of polyhedra the mean volume function is not Morse-Bott in general, despite it being 

522 constructed from the tetrahedral volume. 

523 9 Dynamics 

524 It is easy to check that regular polyhedra are optimal for the gradient of the mean volume 

525 measure. However, not every optimal polyhedron is regular. The results in Section 7 describe 

526 optimal singularities for some classes of polyhedra. By Theorem 4.2 the gradient flow lines 

527 begin and end at singularities. Therefore, if there are only two optimal polyhedra, then the 

528 one with positive mean volume is a sink, while the other one is a source. In order to have 

529 a complete picture about the dynamics of the gradient flow, we would need to address the 

530 singularities with vanishing mean volume and their stable and unstable sets. In particular, 

531 if the stable set of the level singularities are zero-dimensional, then, generically, the mean 

532 volume flow will regularize polyhedra and polyhedral meshes. Unfortunately, this is not so 

533 easy to confirm in general, even though the mean volume is symmetric under orientation 

534 change, and we have an involution N ^ N restricting to an involution of the level set, the 

535 level singularities and the stable and unstable manifolds, respectively. 

536 Let us demonstrate the computation of the level singularities for the pyramid. Consider 

537 '''{^p) = in the case of the pyramid, where p is not coUinear. If 2, 3, 4) = then p2 = Pi 

538 or pi — p3. Suppose that p2 = Pi and v{b, 1, 3) 7^ 0, then by = 1, 3) + i/(5, 1, 4, 3) = 

539 2z/(5, 1, 3) + 4, 3), so that p2 = Pi lies in the plane spanned by pi, p^ and ps, such that 

540 the height of the triangle C(l,4,3) is twice the height of the triangle C(l,5,3). Therefore 

541 the two planar constellations given in Figure 7 are singular pyramids. 




P3 Pi P2 P4 

Figure 7: A family of (planar) singularities for pyramids. 



542 The case j^(1,2,3,4) 7^ is not possible for level singularities. The submanifold of N 

543 where 2, 3, 4) = andpi — p^ is diffeomorphic to the Cartesian product of real projective 

544 3-space with the 3-sphere S0(3) x — RP^ x S^, but it intersects other submanifolds of 

545 the singular set. Therefore the singular set is a union of three submanifolds of N, namely 

546 5^ X and two copies of S0(3) x S^. 

547 In applications, however, we do not start with meshes of negative mean volume. In other 

548 words, if an initial mesh has negative mean volume, then we expect, that there is a problem 

549 with the setup. If we want to apply the flow to polyhedra, they will have well-defined 

550 positive volume, which agrees with the mean volume. If we are given a concrete mesh as a 

551 model of some geometric object, it will approximate the positive volume of the object and 

552 will therefore have positive mean volume. By our definition of regularity, the mean volume 

553 gradient flow regularizes any mesh with positive mean volume, because we avoid the level 

554 set altogether. Even though we cannot rule out the existence of degenerate polyhedra in a 

555 regularized mesh, we have the following useful result. 

556 Theorem 9.1. The mean volume gradient flow regularizes tetrahedra, pyramids and octa- 

557 hedra. The gradient flow of the gradient vector Y given in (7.5) regularizes all hexahedra 

558 p E N satisfying f^ > 0. 
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559 10 Summary 

560 We have shown that the mean volume function or a variation of it can be viewed as a 

561 new, simple and geometrically intuitive regularity and quality criterion, whose gradient 

562 flow improves the quality measured by this criterion. We have proven rigorously, that the 

563 tetrahedron, the hexahedron, the pyramid and the octahedron are regularized by the flow. 

564 Preliminary numerical tests show that the flow of the gradient vector field and its variations 

565 efficiently regularizes the four polyhedron types, which are relevant for the finite element 

566 method. This will be analyzed in detail separately. In the case of the tetrahedron, the 

567 volume function is a Morse-Bott function on S^, whose gradient corresponds to an efficient 

568 smoothing algorithm known as the geometric element transformation method (GETMe), 

569 with one maximal and one minimal 2-sphere, as well as another set of index (2, 2) singular- 

570 ities at level homeomorphic to S*^ x S*^. 

571 

572 Acknowledgements. We thank Joachim Wipper from TWT GmbH Science & Innovation, 

573 Department for Mathematical Research & Services, for fruitful discussions and preliminary 

574 numerical tests. 

575 References 

576 [1] M F Atiyah, R Bott, The Yang-Mills equations over Riemann surfaces^ Philos. Trans. 

577 Roy. Soc. London Ser. A 308 (1983) 523-615 

578 [2] M Atiyah, Collected works. Vol. 5. Gauge theories, Oxford Science Publications, The 

579 Clarendon Press Oxford, University Press, New York (1988) 

580 [3] R Bott, Lectures on Morse theory, old and new, Bull. Amer. Math. Soc. 7 (1982) 

581 331-358 

582 [4] J Cantarella, T Deguchi, C Shonkwiler, Probability Theory of Random Polygons 

583 from the Quatemionic Viewpoint, arXiv:1206.3161 [math.DG] (2012) 

584 [5] R Connelly, I Sabitov, A Walz, The bellows conjecture, Beitrage Algebra Geom. 38 

585 (1997) 1-10 

586 [6] R Connelly, An attack on rigidity. I, II, Bull. Amer. Math. Soc. 81 (1975) 566-569 

587 [7] D E Davies, D J Salmond, Calculation of the volume of a general hexahedron for flow 

588 predictions, AIAA Journal 23 (1985) 954-956 

589 [8] Euclid, Euclid's Elements, Green Lion Press, Santa Fe, NM (2002) 

590 [9] B Griinbaum, Convex polytopes, volume 221 of Graduate Texts in Mathematics, second 

591 edition. Springer- Verlag, New York (2003) 

592 [10] B Kleiner, J Lott, Notes on Perelman's papers, Geom. Topol. 12 (2008) 2587-2855 

593 [11] P M Knupp, Algebraic mesh quality metrics, SIAM J. Sci. Comput. 23 (2001) 193-218 

594 [12] K Kurdyka, T Mostovifski, A Parusinski, Proof of the gradient conjecture of R. 

595 Thom, Ann. of Math. (2) 152 (2000) 763-792 

596 [13] P Mani, Automorphismen von polyedrischen Graphen, Math. Ann. 192 (1971) 279-303 

597 [14] M Morse, The calculus of variations in the large, volume 18 of American Mathematical 

598 Society Colloquium Publications, American Mathematical Society, Providence, RI (1996) 

599 [15] V Rubakov, Classical theory of gauge fields, Princeton University Press, Princeton, 

600 NJ (2002) 

601 [16] IK Sabitov, The volume as a metric invariant of polyhedra. Discrete Comput. Geom. 

602 20 (1998) 405-425 

603 [17] E Steinitz, H Rademacher, Vorlesungen iiber die Theorie der Polyeder unter Ein- 

604 schlufl der Elemente der Topologie, Springer- Verlag, Berlin (1976) 

605 [18] G Strang, G Fix, An analysis of the finite element method, second edition, Wellesley- 

606 Cambridge Press, Wellesley, MA (2008) 



18 



February 26, 2013 



607 [19] D Vartziotis, Analytische Gruppentheorie und Symmetrie-Transformationen, Unpub- 

608 lished manuscript (2008) 

609 [20] D Vartziotis, J Wipper, Fast smoothing of mixed volume meshes based on the ef- 

610 fective geometric element transformation method, Comput. Methods Appl. Mech. Engrg. 

611 201-204 (2012) 65-81 

612 [21] D Vartziotis, J Wipper, M Papadrakakis, Improving mesh quality and finite el- 

613 ement solution accuracy by GETMe smoothing in solving the Poisson equation, Finite 

614 Elcm. Anal. Des. 66 (2013) 36-52 

615 [22] D Vartziotis, J Wipper, B Schwald, The geometric element transformation method 

616 for tetrahedral mesh smoothing , Comput. Methods Appl. Mech. Engrg. 199 (2009) 169- 

617 182 

618 [23] E Witten, Quantum field theory and the Jones polynomial, Comm. Math. Phys. 121 

619 (1989) 351-399 

620 [24] G M Ziegler, Lectures on polytopes, volume 152 of Graduate Texts in Mathematics, 

621 Springer- Verlag, New York (1995) 



19 



February 26, 2013 



